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ABSTRACT 


Continuous  Phase  Frequency  Shift  Keying  (CPFSK)  is  a 
potentially  attractive  modulation  scheme  with  constant 
envelope  and  good  spectral  characteristics  for  energy 
constrained  and  band- limited  satellite  channels. 

This  research  deals  with  Orthogonal  Frequency  Shift 
Keying  (FSK) ,  Minimum  Shift  Keying  (MSK)  which  is  a  special 
case  of  CPFSK,  uncoded  quaternary  CPFSK  and  finally  coded 
quaternary  CPFSK.  Orthogonal  FSK  is  simulated  by  making  the 
modulation  index  (h)  equal  to  one,  and  all  the  other 
simulations  are  performed  with  h=l/2 . 

A  rate  1/2  convolutional  encoder  v/ith  constraint  lengths 
(k)  ,  k=2,  3,  and  4  are  used  in  coded  quaternary  CPFSK 
simulations.  Good  coding  gains  are  obtained  with  only  a  slight 
increase  in  receiver  complexity. 

Soft  decision  with  the  Viterbi  Algorithm  was  applied  to 
all  CPFSK  and  one  MSK  application.  Hard  decision  was  applied 
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I. 


INTRODUCTION 


An  increasing  need  for  communication  between  individuals 
and  nations  is  making  satellite  communication  more  and  more 
attractive  every  day.  Physical  constraints  of  space  based 
applications  force  us  to  choose  bandwidth  and  energy 
efficient  modulation  techniques. 

Continuous  Phase  Modulation  (CPM)  is  a  digital  modulation 
scheme  with  continuous  phase  and  constant  envelope.  Continuous 
phase  makes  it  spectrally  efficient,  and  constant  envelope 
seems  desirable  where  non-linearities  are  present  in  space 
segment  power  amplifiers  and  repeaters. 

Coding  is  used  to  reduce  reception  errors.  The  trade-off 
in  coding  is  a  decrease  in  bandwidth  efficiency.  This  thesis 
mainly  presents  the  research  for  rate  1/2  convolutional 
quaternary  CPM  schemes  with  different  constraint  lengths. 

The  starting  point  for  this  research  was  the  simulation  of 
rate  2/3  trellis  coded  CPFSK.  Tbe  complexity  and  apparent  poor 
performance  of  this  code  forced  us  to  first  study  orthogonal 
FSK,  ana  then  continue  with  MSK,  uncoded  quaternary  CPFSK  and 
finally  finish  with  rate  1/2  trellis  coded  CPFSK.  Due  to  the 
limitation  of  time,  we  were  unable  to  return  to  the  2/3  rate 
codes . 

The  simulatrons  were  done  using  PRO-MATLAB  software  run  on 
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Sun  SPARC  station  computers.  Hardware  limitations  and  the 
complexity  of  the  CPM  schemes  made  this  research  difficult  and 
very  time  consuming. 

Chapter  II  gives  a  brief  review  of  Continuous  Phase 
Modulation .  Chapter  III  focuses  on  Orthogonal  FSK,  and  two 
different  receiver  structures  for  MSK.  Chapter  IV  presents 
uncoded  quaternary  CPFSK  and  coded  quaternary  CPFSK  with 
different  constraint  lengths.  Finally,  Chapter  V  focuses  on 
the  final  conclusions  and  opportunities  for  future  research.. 

An  identical  string  of  fifty  thousand  information  bits  has 
been  used  for  each  simulation  to  provide  a  comparison  betv/een 
systems . 
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II.  BACKGROUND 


This  chapter  focuses  on  the  basic  characteristics  of 
Continuous  Phase  Frequency  Shift  Keying.  In  CPFSK,  frequency 
shifts  are  used  to  transmit  information.  The  phase  of  the 
CPFSK  signal  is  continuous  during  frequency  shifts.  This 
prevents  the  large  spectral  side  lobe  effects  that  are 
observed  in  FSK,  and  PSK.  There  is  also  memory  in  CPFSK  due  to 
the  continuous  phase.  Phase  trellises  will  be  presented  in 
this  chapter  lo.ter  on. 


A.  CPPSK  SIGNAL  DESCRIPTION 

Equation  (2.1)  represents  the  transmitted  signal  in  CPFSK, 
[Ref.  1] 


s  ( t) 


2-^^cos  [27if^t+  0  ( t;  T)  +  ©g] 


(2.1) 


where  Tg  is  the  syiiibul  duration,  the  energy  per  symbol, 

©(t;J)  is  the  time  varying  phase  modulation  induced  by  the 
input  data  sequence,  is  the  carrier  frequency,  and  finally,  0^ 

is  the  initial  phase  of  the  carrier  assumed  to  be  zero  for 
simplicity.  Equation  (2.2)  represents  the  phase  of  the  carrier 
in  the  interval  nT^<.  (n+l)  . 
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(2.2) 


QiC;I)-Kh^  I^->-2Tzhq{  C~nT^) 

K=-«‘ 

The  parameter  n  is  called  the  modulation  index.  Typically  h  is 

chosen  to  be  a  rational  number  h  =  ~  where  m  and  p  are 

P 

relatively  prime  positive  numbers,  in  order  that  takes 

on  a  finite  number  of  states.  The  q(t)  is  a  phase  shaping 
pulse  signal,  I  is  the  input  signal.  For  M-ary  CPFSK,  M=Z'^ 

and  I,,  is  one  of  the  {+1,  +3, . /±{M-1)}  information 

symbols . 

Phase  shaping  signals  can  be  in  a  rectangular  or  raised 
cosine  form.  A  rectangular  form  is  used  in  this  research. 
Phase  trajectories  for  rectangular  pulse  shapes  are  piecewise 
linear.  There  are  p  phase  states  when  m  is  even  and  Pp  phase 
states  when  ni  is  odd.  Smoother  phase  trajectories  can  be 
obtained  by  using  raised  cosine  pulse  shapes. Phase 
trajectories  for  binary  CPFSK  are  illustrated  in  Fig.  2.1. 

B .  ERROR  PERFORMANCE 

Minimum  Shift  Keying  (MSK)  was  used  as  a  baseline 
modulation  scheme;  a  special  case  of  binary  CPFSK  v/ith  h=l/2. 
In  coded  CPFSK  cases  we  also  made  comparisons  with  MSK, 
because,  although  we  send  two  bits  at  a  time,  one  of  them  is 
an  information  bit,  and  the  other  is  a  redundant  cod.'.ng  bit. 
Equation  (2.3)  represents  the  probabiliry  of  error  in  the 
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-2*h*pi 
-3*h*pi 

Fig.  2.1;  Phase  trellis  for  binary  CPFSX. 


presence  of  additive  white  Gaussian  noise  for  MSK.  [Ref.  2] 


Pe  =  Oi 


2^) 


N  "c 


(2.3) 


The  error  performance  can  be  improved  by  increasing  the 
signal  energy,  but  in  satellite  applications  there  are 
limitations  on  power  sources. 

Error  performances  for  each  modulation  scheme  are  given  at 


the  end  of  each  section. 


C .  BANDWIDTH 


Bandwidth 
percent  of  the 
increase  the 
probability  of 


is  the  range  of  frequencies  which  contains  99 
total  signal  power.  A  smaller  value  of  h  will 
bandwidth  efficiency  with  an  increase  in 
error . 
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III.  FSK-MSK 


This  chapter  focuses  on  orthogonal  FSK  and  MSK.  The  Sc'.iae 
modulator  and  bit  mapper  was  used  in  all  simulations.  Input 
bits  were  mapped  to  modulation  symbols  according  to  Table  3.1. 

TABLE  3.1:  MAPPING  FOR  ORTHOGONAL  FSK  AND  MSK. 

input  bit  mapper  (I) 

0 . . -1 

1 . . >  1 


A.  ORTHOGONAL  FSK 

Orthogonal  FSK  simulation  is  accomplished  by  using  binary 
CPFSK  with  h=l.  Coherent  detection  is  used  in  the  demodulator 
and  hard  decisions  are  made  to  detect  the  received  signal.  By 
making  h~l ,  we  have  two  phase  states;  0  and  7t  .  The  phase 
trellis  for  this  system  is  illustrated  in  Fig  3.1.  The  odd 
symbols  represent  the  system  at  state  0,  and  even  s^ndools 
represent  the  system  at  state  7t  .  Multiplying  the  received 

signal  by  offsets  the  total  phase  to  0  at  the  end  of 

each  symbol  and  decreases  the  receiver  complexity.  The 
maximum  value  of  the  correlator  outputs  for  the  two  modulation 
symbols  is  selected  as  the  decision  variable. 
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The  MPiTLAB  source  code  used  for  this  system  is  listed  in 
Appendix  A.  Equation  (3.1)  represents  the  probability  of  bit 


error  for  orthogonal  FSK.  [Ref.  2] 


N, 


(i  .1) 


The  performance  of  this  system  is  illustrated  in  Fig.  3.2.  The 
orthogonal  FSK  simulation  was  done  for  simulation 
verification,  and  as  guidance  to  MSK  simulation. 


Pig.  3.2;  Performaace  of  Orthogonal  FSK 

B.  MSK  WITH  VITERBI  ALGORITHM  DECODING 

MSK  is  a  special  case  of  binary  CPFSK  with  h=l/2  where  the 
phase  changes  ±n/2  at  every  symbol  (see  Fig.  3.3)  .  The  total 
phase  is  ±7t/2  for  odd  symbols  and  the  total  phase  is  0  or  ti: 
for  even  symbols.  When  the  received  signal  is  multiplied  by 
^j2nt/4T,  ^  it  offsets  the  phase  such  that  the  total  phase 

changes  either  by  0  or  by  n  at  every  symbol.  This  reduces  the 
number  of  states  required  in  the  Viterbi  decoder.  [Ref.  3] 
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Fig.  3.3:  Phase  Trellis  for  MSK 

The  state  trellis  for  the  MSK  waveform  is  illustrated  in 
Fig.  3.4.  There  are  two  states,  0  and  n.  There  are  two  paths 
out  of  each  node  corresponding  to  the  two  possible  values  of 
the  input.  Information  bit  1  produces  a  change  in  state  and 
information  bit  0  produces  no  change.  The  Euclidean  distances 
of  the  transmit  symbol  waveforms  from  the  received  signal  are 
used  as  input  to  the  Viterbi  algorithm;  the  Viterbi  algorithm 
chooses  the  path  with  the  minimum  total  distance. 


LO 


Pig.  3.4:  .<3tate  Trellis  for  MSK 

The  cransition  matrix  of  the  MSK  with  Viterbi  algorithm 
decoding  is  illustrated  in  Table  3.2. 

TABLE  3.2:  VITERBI  TRANSITION  MATRIX. 

Ill  0  2  2  1  3 
|l  1  1  2  0  4 

Each  row  of  the  matrix  corresponds  to  a  state.  The  first 
column  of  the  matrix  represent  the  states  that  we  are  coming 
from,  the  second  column  represents  the  input  value,  and  the 
third  column  represents  the  corresponding  number  of  the  point 
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KI)I5 1C  <)<>Jd 


in  the  signal  constellation.  Columns  4,  5,  6  are  the  same 
respectively.  The  MATLAB  source  code  for  this  system  is  given 
in  Appendix  B . 


Theory  vs.  simulation  of  the  MSK  bit  error  probability 
with  Viterbi  detection  is  illustrated  in  Fig.  3.5. 


Error  P  ro  CjO  / 1 1  ty  for  MSK 


Fig.  3,5;  ferf onoance  ot  MSK  with  Viterbi  Detection. 

C.  MSK  WITH  DETECTION  USING  TWO  SYMBOLS  AT  A  TIME 

Equation  (2.1)  can  be  rewritten  as  (3.2) .  [Ref.  4] 


l2 


s  (  t 


(3.2) 


I  Eh 

2  -;;^COS  [  0  (  t)  ]  COS 

'  h 


■2-^sin[©(c)]sin(2Ti;f^t) 

N 


Equation  (3.3)  shows  that  the  phase  of  an  MSK  signal  increases 
or  decreases  linearly  within  the  bit  duration.  [Ref.  4j 

e(t)  =0(O)Tti-;^  (3.3) 

2 

The  minus  sign  in  (3.3)  represencs  information  signal  0,  and 
the  plus  sign  represents  one.  The  signal  phase  can  have  values 
of  ±7r/2  at  odd  multiple  of  or  0  and  it  at  even  multiple  of 

(see  (3.3)  and  Fig.  3.3). 

When  (3.3)  is  inserted  in  (3.2),  it  is  clear  that  the 
input  sequence  during  (0,r^)  affects  only  the  quadrature 

component  of  (3.2),  and  not  the  in-phase  component  over  the 
signalling  interval  (0,2r^)  .  Similarly  the  input  sequence 

daring  (Tj,,2Tb)  affects  only  the  in-phase  component  of  (3.2), 

and  not  the  quadrature  component  over  the  signalling  interval 

(T,3r^)  . 

Table  3.3  was  constructed  from  the  information  discussed 
above.  [Ref.  4] 
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TABLE  3.3:  EFFECT  OF  INFORMATION  BITS  ON  PHASE. 


INFORMATION  BITS 

0(0) 

@(2l) 

0 

n 

n  /2 

0 

0 

-TZ/Z 

1 

0 

nlZ 

1 

TZ 

-TZ  /2 

The  hard  decision  in  the  demodulator  is  done  according  to 
Table  3.3,  and  the  MATLAB  source  code  for  this  system  is 
illustrated  in  Appendix  C. 

Theory  vs.  simulation  of  MSK  bit  error  probability  for 
detection  using  two  symbols  at  a  time  is  illustrated  in  Fig. 

3.6. 


Pig.  3.6;  Performance  of  MSK  for  detection  using  two 
symbols  at  a  time. 
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IV.  CODED  CPFSK  SCHEMES 


Channel  coding  can  be  used  Co  reduce  the  number  of  errors 
when  a  message  is  transmitted  over  a  noisy  channel.  To  achieve 
this  goal,  redundancy  is  added  to  the  information  sequence. 

The  encoder  used  in  this  research  is  a  linear,  sequential, 
finite  state  machine  which  is  implemented  by  using  shift 
registers.  Rate  1/2  encoders  with  constraint  lengths  of  2,  3, 
and  4  are  used.  The  transmitter  block  diagram  for  these 
systems  is  illustrated  in  Fig  4.1,  and  receiver  block  diagram 
is  illustrated  in  Fig.  4.2. 


Fig.  4.1:  Transmitter  Block  Diagram. 
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Pig,  4.2:  Receiver  block  diagram. 

The  bit  mapping  used  in  this  chapter  is  illustrated  in 
Table  4.1.  Two  coded  bits  are  sent  in  each  symbol  duration. 

TABLE  4.1:  MAPPING  FOR  CODED  4-CPFSK. 

input  dibits  mapper  (I) 


00  >  -3 

01  - >  -1 

10  - >  1 

11  >  3 


The  phase  trellis  for  the  uncoded  quaternary  CPFSK  with 
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h=l/2  is  illustrated  in  Fig.  4.3.  By  combining  coding 
techniques  presented  in  this  chapter,  the  probability  of 
choosing  the  wrong  information  sequence  will  be  reduced 


significantly. 


Pig.  4.3:  Phase  trellis  for  quaternary  CPFSK. 

The  demodulator  used  in  this  chapter  is  illustrated  in  Fig 
4.4.  As  is  seen,  the  received  signal  is  first  offset 
and  then  passed  through  four  correlators.  Equation  (4.1) 
represents  the  correlator  outputs. 
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(4.1) 


Pig.  4.4;  Demodulator, 

Real  parts  of  the  correlator  outputs  are  taken,  and  are  used 
as  elements  of  a  four  dimensional  received  signal  V'  :tor.  The 
euclidean  distances  of  this  vector  are  calculated  in  the 
comparator  from  each  of  the  vectOiS  in  the  signal 
constellation.  The  vectors  of  the  signal  constellation  are 
determined  as  follows. 

In  the  absence  of  noise,  the  received  signal  will  be 
e  where  the  phase  state  0^=(O,7r)  and  the 

transmitted  information  sequence  J;^=  ( -3 , -1 , 1 , 3 )  .  Thus  there 
are  eight  vectors  in  the  signal  constellation  determined  by; 
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*r  *  (  t)  s„(  t)  dt] 

(4 .2a) 

(4.2b) 

tj 

“  5 

0 

(4 . 2c) 

r  =¥ 

^kn  *  kn 

(4. 2d) 

Table  4.2  illustrates  the  in  a  tabulated  form.  Each  row 

of  Table  4.2  represents  a  single  point  in  the  signal 
constallaticn . 


TABLE  4.2:  SIGNAL  CONSTELLATION. 


^k/^n 

-3 

-1 

1 

3 

-3  (0^=0) 

1 

0 

0 

0 

-3  (e^=n) 

'1 

0 

0 

0 

-1  (0^=0) 

0 

1 

0 

0 

-1  (0(.=-rt) 

0 

-1 

0 

0 

1  (©j^-O) 

0 

0 

1 

0 

1  (0jt  =  Tt) 

0 

0 

-1 

0 

3  (0^=0) 

0 

0 

0 

1 

3  {e^=n) 

0 

0 

0 

-1 

Equation  (4.1)  is  represented  as  (4.3)  in  the  MATLAB 
program. 

C^=  x’(k\t)  s„{k^t)  (4.3) 

A-O 


19 


where  At  is  the  sampling  period,  and  T^/Ac  is  the  number  of 


samples  per  symbol  (N)  .  In  MA.TLAB  simulation,  Table  4.2 
becomes  Table  4.3. 


TABLE  4.3 j  SIGNAL  CONSTELLATION  FOR  MATLAB  APPLICATION. 


In 

r  .  ■  ' 

-3 

- 1 

1 

3 

-3 

o 

11 

© 

N 

0 

1 

0 

-3 

II 

© 

-N 

0 

-1 

0 

- 1 

(0;,  =  O) 

0 

N 

0 

1 

- 1 

© 

II 

0 

-N 

c 

-1 

1 

1 

(0^=0) 

1 

0 

N 

0 

1 

(0;.-'n:) 

-1 

0 

-N 

0 

3 

(0.^=0) 

0 

1 

0 

N 

3 

© 

II 

0 

-1 

0 

-N 

A.  CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  TWO 

The  coding  used  in  this  section  is  created  by  using  two 
shift-registers.  [Raf.  5]  The  encoder  is  illustrated  in  Fig. 
4.5,  and  the  state  trellis  for  this  system  is  illustrated  in 
Fig  4.6  (b) . 

The  demodulator  calculates  the  euclidean  distances  of  the 
received  signal  vector  from  the  eight  vectors  illustrated  in 
Table  4.4.  At  each  symbol  duration,  eight  distance  values  are 
entered  to  the  Viterbi  algorithm  and  the  true  path  is  chosen 
according  to  minimum  distance. 


2  0 


Most  significant  bit 


Least  significant  bit 


Fig,  4.5:  Encoder. 


States 


Fig.  4.6:  Phase  Trellis  (a),  and  State  Trellis  for  coded  4 
CPFSK  with  k*=2  (b)  . 


When  comparing  Fig.  4.3  and  Fig  4.6  (a),  it  is  clear  that 
at  each  symbol  duration,  the  probability  of  choosing  the  wrong 
state  is  reduced  by  this  simple  coding.  Later,  in  the 
following  sections,  more  complicated  codes  will  be  studied  and 
the  coding  gains  that  are  obtained  will  be  seen  more  clearly. 

TABLE  4.4:  SIGNAL  CONSTELLATION. 


o 

r-l 

1 

O 

0,1,0,17 

1 ,  C ,  17,  0 

-1,0, 

~N,  0 

|o,  -N,0,  -1 

0 , 17,  0 , 1 

A/,0,1,0 

-17,  0, 

-1,0 

The  transition  matrix  for  this  system  is  illustrated  in  Table 

4.5. 


TABLE  4.5;  VITERBI  TRANSITION  MATRIX. 


1 

0  1 

1 

0 

2 

2 

1 

3 

2 

1 

4 

1 

2  5 

1 

2 

6 

2 

3 

7 

o 

3 

8 

Each  row  of  the  matrix  corresponds  to  a  state.  The  first 
column  of  the  matrix  represent  the  states  that  we  are  coming 
from,  the  second  column  represents  the  input  value,  and  the 
third  column  represents  the  corresponding  point  in  the  signal 
constellation  that  is  illustrated  in  Table  4.4;  the  points  are 
numbered  rowwise.  All  other  columns  are  the  same  respectively. 

The  performance  of  this  system  is  illustrated  in  Fig.  4.7. 
The  coding  gain  obtained  will  not  be  worth  what  must  be 
sacrificed  in  bandwidth  efficiency  (see  Fig.  4.7)  .  Since  this 
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was  the  initial  experiment  in  the  use  of  coding,  it  provided 
a  baseline  for  cur  study.  The  MA.TLAB  source  code  for  this 
system  is  listed  in  Appendix  D. 


E-ft-o"  r-Q  t>o  Di  tity  J'or*  MSK  0*^0  ^ — OPF’SK  k^2 


Fig  4.7:  Performance  of  the  coded  4-CPFSK  with  ]t=2  . 

B,  CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  THREE. 

The  coding  used  in  this  section  was  created  by  using  three 
shift- registers .  The  encoder  is  illustrated  in  Pig.  4.8,  and 
the  state  trelli.s  for  this  system  is  illustrated  in  Fig  4.9. 
[Ref.  5] 

The  demodulator  calculates  the  euclidean  distances  of  the 
received  signal  vector  from  the  eight  vectors  presented  in 
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Pig  4. 


rig  A. 


8  s  Encoder , 


States 


9  State  Trellis  for  coded  4-CPFSK  with  k=3 „ 
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Table  4.4.  At  each  symbol  duration,  eight  distance  values  are 
entered  to  the  Viterbi  decoder  and  the  true  path  is  chosen 
accc ■' ding  to  minimum  distance. 


The  transition  matrix  for  this  system  is  given  in  Table 


4.6. 

TABLE  4.6s  VITERBI  TRANSITION  MATRIX. 

101102225226 
33  7  338413414 
202201126125 
314313438437 


The  performance  of  this  system  is  illustrated  in  Fig. 
4.10. 


Error  «=»roboe>'«*ty  ^or  MSK  ori<i  4,— CP»rSK  with  U  3 


Fig.  4. 1C:  Perfonrance  of  coded  4-CPFSK  with  k=3 . 
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The  coding  starts  working  efficiently  after  6 . 7  dB 

At  this  signal  -  to-noise  ratio  almost  1  dB  gain  is  obtained  and 
at  7.6  dB  e:^/Nq,  the  coding  gain  reached  almost  2  dB. 

Signal  -  co-noise  ratios  for  8.8  dB  and  10  dB  were  also 
tested  and  no  errors  were  detected  in  fifty  thousand 
information  bits. 

The  MATLAB  source  code  for  this  system  is  listed  in 
Appendix  E . 

C.  CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  POUR. 

The  final  coding  studied,  uses  four  shift  registers.  The 
encoder  is  illustrated  in  Fig.  4.11,  and  the  state  trellis  for 
this  system  is  illustrated  in  Fig.  4.12.  [Ref.  5] 

In  this  application,  the  demodulator  calculates  the 
euclidean  distances  of  the  received  signal  vector  from  eight 
different  vectors  presented  in  Table  4.4.  At  each  symbol 
duration  eight  distances  are  entered  to  the  Viterbi  decoder, 
and  the  true  path  is  chosen  according  to  minimum  distance. 

The  transition  matrix  of  this  system  is  illustrated  in 
Table  4,7  and  the  MATLAB  source  code  is  listed  in  Appendix  F. 

The  performance  of  this  coding  is  illustrated  in  Fig. 
4.13.  The  coding  starts  working  efficiently  at  5.9  dB  E^^/N^. 

At  this  signal  - to-noise  ratio  almost  l  dB  gain  is  obtained  and 
at  6.8  dB  Ej^/Eq,  the  coding  gain  reached  almost  2  dB. 
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4.11:  Encoder. 

1 

2 

3 

4 

5 

6 

7 

8 

4.12:  Trellis  for 
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coded  4-CPFSK  with  k=4. 
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TABLE  4.7:  VITERBI  TRANSITION  MATRIX 


1 

0 

1 

n 
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Signal - to-noise  ratios  of  7.7  dB,  8.8  dB,  and  10  dB  are 
also  tested  and  no  errors  were  detected  for  fifty  thousand 
information  bits. 


CrrQr  ProooOMIty  for  MSK  ono  vvithv  kw  4 


Fig.  4,13;  Performance  of  the  coded  4-CPPSK  with  ks=4. 


28 


V,  CONCLUSIONS  AND  FURTHER  RESEARCH 


This  thesis  mainly  presented  the  simulation  results  for 
rate  1/2  convolutional  quaternary  CPM  schemes  with  different 
constraint  lengths.  Using  different  codes,  the  minimum  merge 
length  increased,  and  larger  minimum  Euclidean  distances  were 
obtained.  Better  error  performances  were  obtained  as  a  result 
of  these  codes. 

This  research  can  be  extended  in  several  ways.  First,  rate 
1/2  codas  can  be  improved  by  using  different  constraint 
lengths,  and  shift  register  connections.  The  same  encoders  can 
also  be  tested  with  different  modulation  indexes. 

More  complicated  rate  1/2  codes  with  h=l/2,  can  be 
simulated  very  easily  by  making  small  changes  to  the  MATLA3 
source  codes  that  are  listed  in  Appendices  D  thru  F.  The 
signal  constellation  that  is  presented  in  Table  4.4  can  be 
used  for  every  new  code. 

Second,  rate  2/3  codes  tor  different  constraint  lengths, 
shift  register  connections,  and  modulation  indexes  can  be 
studied,  and  best  systems  with  minimum  receiver  complexity, 
and  maximum  power -bandwidth  performances  can  be  found. 
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APPENDIX  A 


%Thi3  program  creates  a  m  by  n  binary  input  matrix. 

function  [random_matrix] =random(m, n) ; 

a=rand (m, n) ; 

b=ones (m, n) *0 . 5 ; 

random  matrix=f loor (a+b) ; 


%  This  program  d,s  used  as  mapper.  It  converts  0  to  -1  and  1 
%  to  0 . 

function  [mapper_output] =mapper ( input ) ; 

[m,a] =si2e (input) ; 
for  x=l:a 

y=ones (m, 1) *2^ (2 -x) ; 
z= [z,y] ; 
end; 

mapper_output= input . *z ; 
mapper_output=mapper_output ' ; 
mapper_output=  (mapper__output)  -1; 


%  This  program  is  used  as  CPFSK  modulator.  Inputs  to  this 
%  program  are  "T=syrtibol  duration",  "h=modulation  index",  and 
%  "modulator_input=mapper  output". 

function  [mod__output]  =modulator  (T,  h, modulator_input )  ; 

z=-length  (modulator_input)  ; 

teta0=0 ; 

for  sample=l:z 

I=modulator__input  (sample)  ; 
if  saLmple==i 

teta=0;  %Sets  the  initial  phase  to  zero, 
else 

tetaO=tetaO+modulator_input (sample- 1) ; 
teta=pi*h*tetaO ; 
end 

time=0 : 0.0005 :T' 0.0005, • 
time=time+ ( sample- 1) *T; 
xx=pi*I*h* (time/T- (sample- 1) ) +teta; 
mod_output (sample, : )= (exp(i*xx) ) ; 

end 
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%  This  program  is  used  to  add  Gaussian  noise  to  modulator 
%  output,  Inputs  to  this  program  are  "m  as  modulator  output", 
%  "coeff  as  standard  deviation  of  the  noise",  and  "T  as  symbol 
%  duration".  The  constant  0.005  as  sampling  period. 

function  [awgn_noise]  =av7gn_ch  (m,  coeff ,  T) 
rand ( ' seed' ,99); 
rand ( ' normal ' )  ; 

awgn__noise=coef  f  *  ( rand  {m,T/ 0.005)  +  j  *rand  im,  T/0 . 005 )  )  ; 


%  This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "mod_output  as  received  signal".  In  this  demodulator  hard 
%  decision  is  used. 

function  [demod_output] =demodulator (T, h, mod_output) ; 
map= [ - 1  1] ; 

[m,  a] =size (mod_output ) ; 

Cime=0 :0.0005:T-0.0005; 
for  sample=i:m 
x=rem (sample, 2 ) ; 
phase__coef  f  =  [  (l-x)  *pi]  ; 
phase=exp (i*phase_coef f ) ; 
for  m_ary=i:2 

xx=exp { i*time*pi*h*map (m_ary) /T) ; 
match=xx*conj (mod_output (sample,  : ) ' ) ; 
match=phase*match; 
match_ouc= [match_out  match] ; 

end 

[y, i] =max(real (match_out) ) ; 
demod_output ( sample ) =map ( ceil ( i/1 )  )  ; 
macch_out= [  ]  ; 
end 
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function  [Nt),I,e]  =check{x,y) 

%  SYMBOL  ERROR  CHECK 

%  Paul  H.  Moose 

%  Naval  Postgraduate  School 

%  09-01-91 

o, 

o 

%  This  m-file  locates  the  positions  in  vectors  x  and  y  that  do 
not 

%  agree.  It  returns  a  one  in  e  if  they  do  not  agree  and  a  zero 
if  they 

%  do  agree,  e  is  the  error  vector  if  x  and  y  are  binary.  I  is 
a  vector 

%  of  error  location  numbers.  Nb  is  the  sum  of  the  elements  of 
e . 

e= (x~=y) ; 

I=f ind (e) ; 

Nb=sum(e) ; 


This  program  is  used  as  master  program.  It  calls  different 
functions  and  simulates  whole  communication  system. 


T=0.015 

h=l; 

m=50000 

n=l; 


(I  (p  II 
II  h" 

"m" 


is  used  as  symbol  duration, 
is  used  as  modulation  index. 


and  "n"  are  used  to  create  binary  matrix. 
coeff=0;  %  "coeff"  is  used  as  standard  deviation  of  the  noise. 
[input_matrix] =random(m,n) ; 

[mapper_output] =mapper ( input_matrix) ; 

[m.od_output]  =modulato  (T,  h,mapper_output)  ; 

[awgn_noise]  =awgn__ch  (m,  coeff ,  T)  ; 
mod_output=mod_output+awgn_noise; 

[dsmcdulator_output j  -demodula  (T,  h,rriod_outpuL)  ; 
error=check (mapper_output , demodulator_output ) ; 
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APPENDIX  B 


%  This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "mod_outpat  as  received  signal".  In  this  demodulator  hard 
%  decision  is  used. 

function  [demod_output] ^demodulator (T, h,mod_output ) ; 
map=  I  - 1  1 ]  ; 

Im, a] =size (mod_output ) ; 
time=0 : 0 . 0005 ; T; 
for  sample=i:m 

mod=mod__outpuc  (sample,  :)  *exp  (j  *(  sample- 1)  *pi/2 )  ;  %Offsets 

%the  phase 

for  m_ary=l:2 

xx=exp ( i*time*pi*h*map (m_ary ) /T) ; 
match=xx*conj  (mod_output  (saunple,  : )  '  )  ; 
inatch__out=  [match_out  match]  ; 

end 

demod_output  (sample,  : )  ==  (real  (match_out)  )  ; 

match_out=  [  ]  ; 

end 

R-  ( j  ♦demod__output  ( : ,  1 )  +demod_output  (;,2))/31; 

[demod_output]  =eucdis  (2 , R)  r  %Calculates  the  euclidean  distance 

%of  the  received  signal 


function  D  =  eucdi3(q,R) 

%  EUCLIDEAN  DISTANCE  METRICS 

%  Paul  H.  Moose 

%  Naval  Postgraduate  School 

o-  c\iZ  n  no 

% 

%  This  M-file  finds  Euclidean  distance  of  elements 
%  in  vector  R  from  2'^q  unit  amplitude  vectors 

equally  spaced  on  the  unit  circle.  It  stores  these  as  rows 
of  D . 

N=2^q; 

L=length (R) ; 
index=l :N* 
dph=2*pi/N; 

M0=exp ( j  *dph. * ( index- 1) ) ; 
for  1=1:L 

D  ( 1 ,  -.  )  =abs  (R  (1)  .  *ones  (MO)  -MO)  ; 
end 
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function.  PHN  ==  sof  tv  (k,  K,  Np,  PH,  T,  U) 

%  Soft  Viterbi  Decoder 

%  Paul  H.  Moose 

%  Univ.  degli  Studi  di  Padova 

%  17-05-91 

a 

"o 

%  This  M-file  decodes  k  bit  msgwords  from  2'^n  real  metrics 
%  (These  may,  for  exar,  represent  the  "distance"  of  the 
%  received  modulation  value  from  each  of  2^n  modulation 
values . ) 

%  The  state  transition  information  for  a  2'^K  state  trellis  is 
in 

%  the  2'^K  by  3*2^k  matrix  T.  Each  of  the  2'^k  entering  paths 
to 

%  each  state  has  its  source  state  (one  of  2''K)  ,  path  msgword 
( one 

%  of  2^k)  and  path  codeword  (one  of  2^n)  listed  in  the  state 
row. 

%  The  path  histories  are  kept  in  matrix  PH  that  is  2'^K  by 
3*Np. 

%  The  path  history  for  each  state  contains  source  state,  path 
%  weight  and  path  codeword  for  Np  previous  states. 

%  The  output  PHN  is  the  update  of  PH,  the  new  path  history. 
%  The  decoded  code  word  is  in  the  last  column  of  PHN.  (They 
should 

V  Ul^X  • 

%  The  past  histories  are  undated  on  the  basis  of  the 
"minimum 

%  metric".  You  can  change  this  to  the  "rna,ximum  metric"  if 
desired  as 

%  indicated  in  the  comments  in  the  code. 

% 

for  j=l:2'"K 

X( j ,2) =D (T( j , 3) ) +PH{T( j , 1) ,2) ;  %path  weight 
X  ( j  ,  1)  =T  ( j  ,  1)  ;  %path  SOI’ roe  state 

X  ( j  ,  3  )  =T  { j  ,  2)  ;  %path  code  word  T  ( j  ,  3 )  ,  Chg  to  T  ( j  ,  2 )  for 

insf^word . 

for  l=2:2"k 

wt  =  D(T(j,3*l))  +PH(T(_i  ,3*1-2)  ,2)  ; 
if  wt  <  X(j,2)  %The  <  selects  min  metric 
X(j  ,2)  =wt; 

X{j,l)=T(j, 3*1-2)  ; 

X(j, 3)=T(j, 3*1-1)  ;  %Chg  to  T(j,3*l)  for 

codeword. 

else 

end 

end 

%  VJe  now  need  to  append  old  paths  to  new  paths  to  get 
survivors . 

PHN(j,:)=  [X(j,;)  PH  (X  ( j  ,  1 )  ,  1 :  3*Np- 3  )  ]  ; 
end 
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%  This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  Scirae  as  in  Appendix  A. 

clear 

temp= [1  0121  4 ; 1  203];  %Is  used  as  Viterbi  pat  i 

%matrix. 

T=0 .015; 
h=l/2; 
m=50000 ; 
n=l; 

coef f =0 ; 

[input_matrix] =random(m,n) ; 

[mapper_output] =mapper (input_matrix) ; 

[mad_output] -modulate (T, h,mapper_output) ; 

[awgn_noise] =awgn_ch(m, coeff ,T) ; 
mod_ou  t  pu  t  =inod_ou  t  pu  t + awgn_no  i  s  e  ; 

[demodulator__output]  =demodula  (T,  h,  mod_output )  ; 
viterbi_path__matrix=temp; 

TT=zeros (2,60) ; 
for  x=l:m 

D=demodulator_output (x, : ) ; 

[TT] =3oftv(l, l,20,TT, viterbi_path_matrix,D) ; 
received_signal (x)=TT(i,60) ; 
end 

input_matrix=input_matrix(l9 :m) ; 
received_signal*received_signal (l:m-19) ; 
error=check(received__signal,  input_matrix)  ; 
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APPENDIX  C 


This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
"in  as  received  signal".  In  this  demodulator  hard  decision 
is  used. 

function  [demodulator^output] =demodula (T, h, in) 

N=T/.0005; 

[k,  a]  =size  (in)  ; 

P  =l:a; 
y3=-l; 
for  m=3 : k 

phil  =  exp ( j  *pi* (m- 1) /2 ) *exp ( j  *  (pi/2)  . * (p/N) ) ; 
phi2=exp ( j*pi* (m-2) /2) *exp ( j* (pi/2) . * (p/N) ) ; 
x=rem (m, 2 ) ; 
if  (x==0) 

y=sum(real(in(m-2, :) ) .*abs(imag(phi2) ) )  + 

sum (real (in(m-l, : ) ) . *abs (imag(phii) ) ) ; 
if  (  (y<0  &  y3>0)  I  (y>0  &  y3<0) ) 
d=l; 

else 

d=-l; 

end 

else 

y  =  sum(imag(in(m-2,  :)  )  .*abs(real  (phi2)  )  )  +• 

sum ( imag ( in (m- 1 , ; ) ) . *abs ( real (phil ) ) ) ; 
if  (  (y>0  &  y3<0)  I  (y<0  &  y3>0) ) 
d=  - 1  ; 

else 

d=l; 

end 
end 
y3=y ; 

demodulator_output (m) =d; 
end 
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%  This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  same  as  in  Appendix  A. 

clear 
T=0 . 015; 
h=l/2; 
m-  5  0  0  0  0 ; 
n=l; 

coef f =0 ; 

[input_matrix] =random(m,n) ; 

[mapper_output] =mapper (input_matrix) ; 

[mod_output]  =iTiodulato  (T,h,mapper_output)  ; 

[awgn_noise] =awgn_ch (m, coef f ,T) ; 
mod_ou  tpu t  =mod_ou  t  pu  t  +awgn_no i s  e ; 

[demodulator__output]  =demodula (T, h,mod_output)  ; 
error=check;  (mapper_output  (2  :m-2)  ,  demodulator_output  (4  ;m)  )  ; 
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APPENDIX  D 


%  This  program  axmulates  the  encoder  that  is  given  in 
%  Fig.  4.4.  The  input  to  this  program  is  a  m  by  1  binary 
%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 
%  output  .  First,  two  shift  registers  are  initialized  and  then 
%  input  are  shifted. 

function  [ tcm_output ] =tcm { input_matrix) ; 

[m, a] =size ( input_macrix) ; 

f  f  1=  [0;  input_ma.trix  ( :  ,  1)  ]  ;  %  Intializes  the  first  flip-flop; 
f f l=f f 1 { 1 :m, 1 ) ; 

f f2= [0 ; f f 1]  ;  %  Intializes  the  second  flip-flop; 

.Cf2  =  ff2  (l:m,  1)  ; 

tcm_output { : , 1) =f f 1 ;  %  Most  significant  bit. 

tcm__output  ( : ,  2 )  =f  f  2  ;  %  Least  signif  icicent  bit. 


%  This  program  is  used  as  mapper.  It  simulates  the 
%  Table  4.1. 

function  rinapper_output]  =mapper( input)  ; 

[m, a] =size (input) ; 
for  x=l:a 

y=ones  (m,  1)  *2^^  ( 3  -  x)  ; 
z=  [z,y]  ; 
end 

mapper__output  =  input .  *z; 
mapper_output=iaapper_output '  ; 
mapper__output=snm  {mapper_output )  -  3 ; 


%  This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "mod_output  as  received  signal". 

function  [demod_output] =demodulator (T, h,mod_output) ; 
mapl=  [  -  3  - 1  3]; 

[m, a] =size {mod_output ) ; 
time=0 : 0 . 0005 : T; 
for  sample=l:m 

mod=exp { j  * (sample- 1) *pi/2) *mod_output (sample, : ) ; 
for  m_ary=l:4 

xx=exp (i*time*pi*h*mapl (m_ary) /T) ; 
matchl=xx*conj (mod' ) ; 
match_out=  [match__ouc  matchi'  ]  ; 

end 

demod_oucput (sample, : ) =distance (real (match_out) ) ; 

match_out= [] ; 

end 


%  This  program  calculates  the  euclidean  distances  of  the 
%  received  signa.l  from  signal  constellation. 

function  [dist] =distance (R) ; 

matrix=ro  -1  0  -31  %  This  matrix  represents  the  signal 

0  1  0  31  %constellation, 

1  0  31  0 
-1  0  -31  0 
0  -31  0  -1 
0  31  0  1 
31  0  1  0 
-  3 1  0  - 1  0  ]  ; 
d=ones (8,1)  ; 
dQ-d*R; 

dd=  (dd-matrix)  .'^2; 
dist=sum (dd' ) ; 
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%  This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  ocher  functions  are  same  as  in  Appendix  A. 

clear 

temp=[l  01102213214 

12512623  7  23  8]; 

T=0 . 015 ; 
h=l/2; 
m=  5  0  0  0  0 ; 
n  =  l; 

coef f =0 ; 

[ input__matrix]  =random (m,  n)  ; 

[tcm_output] =tcm(input_matrix) ; 

[mapper_output] =mapper ( tcm_output) ; 

[mod_output] =modulato (T, h,map) ; 

[awgn_noise] =awgn_ch (m, coef f , T) ; 
mod_output  =mod_ou tput +awgn_noise ; 

[demodulator_output] =demodula {T,h,mod_output) ; 
viterbi_path_matrix-cemp; 

TT=2eros (2,120) ; 
for  x=l:m 

D=demodulator _output (x, : ) ; 

[TT] =sof tv (2,1,40, TT, vicerbi_path_matrix, D) ; 
received_signai  (x)  =-TT  (1,120)  ; 
end 

ww=mb (2 , received_signal ) ; 
ww=ww  U  :  2  .*  2  *m)  ; 

error=check  ( inpuc_matrix(l  :m-41)  ' ,  ww  (42  ;m) )  ; 
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APPENDIX  E 


%  This  program  simulates  the  encoder  that  is  given  in 
%  Fig.  4. .  The  input  to  this  program  is  a  m  by  1  binaiy 
%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 
%  output  .  First,  the  three  shift  registers  are  initialized 
%  and  then  input  is  shifted. 

function  [tcm_output] =tcm(input_matrix) ; 

[m,a]  =:size  (  input_matrix)  ; 

f f 1= [0 ; input_matrix ( : , 1) ] ;  %  Initialize  the  first  flip-flop, 
f f l=f f 1 

f f 2=  [0 ; f f 1] ;  %  Initialize  the  second  flip-flop. 

ft2=ff2 (l:m, 1) ; 

f f 3= [0; f f2] ;  %  Initialize  the  third  flip-flop, 

f f 3=f f3 (1 :m; 1) ; 

tci,Ti_output  ( :  ,  1)  =abs  (abs  (f  fl-ff2)  -f  f3)  ;  %  Most  significant  bit . 
tcm_output ( : , 2) =f f2;  %  Least  significant  bit. 
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funccd.on  [PHN]  =  viterbi  (k,  Np,  PH,  T,  D) 

%  Soft  Viterbi  Decoder 

%  Paul  H.  Moose 

%  Mercury  Digital  Communications 

%  06-09-91 

*0 

%  This  M-file  decodes  k  bit  msgwords  from  2'^n  real  metrics 
%'  (These  may,  for  example,  represent  the  "distance"  of  the 
%  received  modulation  value  from  each  of  2'^n  modulation 
values . ) 

%  The  state  transition  information  for  a  2^K  state  trellis  is 
in 

%  the  2^K  by  3*2'^k  matrix  T.  Each  of  the  2^k  entering  paths 
to 

%  each  state  has  its  source  state  (one  of  2''K)  ,  path  msgword 
(one 

%  of  2'^k)  and  path  codeword  (one  of  2'^n)  listed  in  the  state 
row. 

%  The  path  histories  are  kept  in  matrix  PH  that  is  2''K  by 
3*Np. 

%  The  path  history  for  each  state  contains  source  state,  path 
%  weight  and  path  codeword  for  Np  previous  states. 

%  The  output  PHN  is  the  update  of  PH,  the  new  path  history. 
%  The  decoded  msg  word  is  in  the  last  column  of  PHN.  (They 
should 
%  "merge". 

%  The  past  histories  are  undated  on  the  basis  of  the 
"minimum 

%  metric".  You  can  change  this  to  the  "maximum  metric"  if 
desired  as 

%  indicated  in  the  comments  in  the  code. 

% 

P=FH(: ,2) ' ; 
wt=T(: ,3:3;3*2"k) ' ; 
ux=T(: ,l:3:3*2^k-2) ' ; 
aa=D (wt ( ; ) ) ; 
bb=P (ux (:)); 

wt ( ; ) =aa+bb;  %This  contains  all  weights (columns )  for  each 
state ( row) 

[a, b]  =min (wt )  ;  %  Use  maJi:(wt)  here  for  maximum 
X(  :  ,  2)  =a'  ; 

X(: , 1) =diag(T( ; ,3.*b-2' ) ; 

X ( : , 3 ) -diag (T ( : , 3 . *b- 1) ) ; %path  msgword.  Chg  to  3 . *b  for 
codeword 

%We  now  need  to  append  old  paths  to  new  paths  to  get 
survivors . 

PHN=[X  PH{X( ;  ,  1)  ,  l:3*Np-3) ]  ; 
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%  This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  same  as  in  Appendix  A. 

clear 

temp=[l  01102225226 
337338413414 
202201126125 
31431343843  7]; 

T=0.015; 
h=l/2; 
m=50000 ; 
n=l; 

coef f =0 ; 

[input_matrix] =random{m,n) ; 

[tcm_output] =tcm(input_matrix) ; 

[mapper_output] =mapper {tcm_output) ; 

[mod_output] =modulato (T, h,map) ; 

[awgn_noise] =awgn_ch (m, coef f ,T) ; 
mod_output=mod_output+awgn_noise; 

[demodulator_output] =demodula (T,h,mod_output) ; 
viterbi_pat  h_ina  t  r  ix=  t  emp  ; 

TT= zeros  (4,120) ; 
for  x=l:m 

D=demodulator_output  (x, : ) ; 

[TT] =viterbi (2,40, TT, viterbi_path_matrix,D) ; 
received_signal (x) =TT(1, 120) ; 
end 

ww=mb  (2,  received__signal)  ; 
ww=ww(l:2  :2*m)  ; 

error=check (input (l:m-41) ' ,ww(42:m) ) ; 


43 


APPENDIX  F 


%  This  program  simulates  the  encoder  that  is  given  in 
%  Fig.  4.  The  input  to  this  program  is  a  m  by  1  binary 
%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 
%  output  .  First,  the  four  shift  registers  are  initialized 
%  and  then  input  is  shifted. 

function  [tcm_output] =tcm(input_matrix) ; 

[m, a] =size ( input_matrix) ; 

f f 1= [0 ; input_matrix ( : , 1)  ]  ; 

ffl=ffl {l:m, 1) ; 

ff2=  [0;ffl]  ; 

ff2=f f2 (l;m, 1) ; 

ff3=  [0;ff2]  ; 

ff3=ff3 (l:m, 1) ; 

ff4=  [0; ff3]  ; 

ff4=ff4 (l:m, 1) ; 

tcm_output ( : , 1) =abs (f f 1-abs (ff2-f f4) ) ; 
tcm_output  h / 2 ) =f f 3 ; 
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%  This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 


%  other  : 

functions 

are 

(  same 

as 

1  in  Appendix  A 

clear 

temp= [ 

1 

0 

1 

1 

0 

2 

2 

2 

6 

2 

2 

5 

4 

3 

7 

4 

3 

8 

3 

1 

3 

3 

1 

4 

5 

2 

5 

5 

2 

6 

6 

0 

2 

6 

0 

1 

7 

3 

7 

7 

3 

8 

8 

1 

4 

8 

1 

5 

1 

2 

5 

1 

2 

6 

2 

0 

2 

2 

0 

1 

3 

3 

8 

3 

3 

7 

4 

1 

3 

4 

1 

4 

5 

0 

2 

5 

0 

1 

6 

2 

5 

6 

2 

6 

7 

1 

4 

7 

1 

3 

8 

3 

7 

8 

3 

8]  ; 

T=0 .015 

/ 

h=l/2; 

m=50000 

/ 

r.=  l; 

coef f =0 

; 

[input_ 

matrix] : 

=random{m,: 

n) 

/ 

[tcm_output] =ccm(input_matrix) ; 

[mapper_output]  =mapper  ( tcm__QUtput)  ; 

[mod_output] =modulato (T,h,map) ; 

[awgn_noise] =awgn_ch (m, coeff ,T) ; 
mod_output=mod_output+ a wgn_no i s e ; 
[demodulator_output] =demodula (T, h, mod_output ) ; 
viterbi_path_matrix=temp; 

TT=zeros (8,120)  ; 
for  x=l:m 

D=demodulacor_output (x, : } ; 

[TT]  =viterbi  (2, 40  ,  TT,  vitert»i_path_matrix,  D)  ; 
received_signal (x) =TT (1, 120) ; 
end 

ww=mb (2 , received_signal) ; 
ww= ww ( 1 ; 2 : m*  2 )  ; 

errar=check ( input_matrix ( 1 ;m-42) ' , ww (43  ;m)  )  ; 
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